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Abstract. Internal friction effects are responsible for line widening of the 
resonance frequencies in mechanical oscillators and result in damped oscillations 
of its eigenmodes with a decay time Q/lu. We study the solutions to the equations 
of motion for the case of spherical oscillators, to be used as next generation 
of acoustic gravitational wave detectors, based on various different assumptions 
about the material's constituent equations. Quality factor dependence on mode 
frequency is determined in each case, and a discussion of its applicability to 
actual gravitational wave detectors is made on the basis of available experimental 
evidence. 



PACS numbers: 04.80.Nn, 95.55.Ym 



1. Introduction 

Spherical Gravitational Wave (GW) detectors will almost certainly be, in one of 
its variants (solid, hollow or dual), the next generation of acoustic antennae, due 
to their multimode-multifrequency capabilities ^, ^, H, ||, 0, as well as 

their potentially enhanced sensitivity relative to their currently operating cylindrical 
counterparts j|, ^, 0, 12 1. Conviction that this is going to be the case has 



encouraged a remarkable research effort within the GW community, and a variety 
of important topics have been addressed, ranging from theoretical to practical aspects 
of the problem. Several countries worldwide are currently developing projects to build 
and operate spherical detectors |l^ . 

Some of the salient properties of a spherical body as a detector of GWs can be 
accurately established by means of an idealised model, where dissipative effects can be 
safely overlooked |16, ij, Such effects are however very important when it comes 



to determining the ultimate sensitivity of the device, due to their fundamental bearing 
on the noise characteristics of the antenna via the fluctuation-dissipation theorem, 
more specifically its spectral density. 

In a solid elastic body there are many oscillation modes in addition to the 
fundamental one chosen for GW detection, which therefore add noise to the frequency 
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band of operation. It has been suggested ^ that normal mode analysis of this 
contribution to the noise level can be either inconvenient or too expensive in terms of 
computational cost if heat is e.g. generated by a laser spot on the oscillating solid, such 
as happens in an interferometric GW detector or in a dual sphere Q. Experimental 
measurements clearly show that this is the case 

The naive approach, and actually the one most often taken, to account for 
dissipative effects is to add an ad hoc term proportional to the velocity in the solid's 
eigenfrequency equation. This results in a damped oscillatory behaviour of the sort 

e-'^*/'5 sinut (1) 

where the quality factor Q is assumed constant (i.e., time independent), and accounts 
for the linewidth of the mode with frequency uj. It is well known that this quality 
factor is different for different oscillation modes, so the set of Q's is very large, and 
the procedure appears highly artificial. On the other hand, there has been criticism 
of the concept of dissipation forces proportional to velocity in terms of actual noise 
spectrum in an elastic solid , and suggestion that a better description may rather 
derive from the use of complex frequencies. 

In the specialised literature on the subject (see and references therein), 
viscoelastic effects are often described by means of so called constituent equations. 
These are extensions of Hooke's law relating the stress and strain tensors in the solid. 
Much like in a simple one-dimensional spring, internal friction forces can be considered 
proportional to the instantaneous velocity of the oscillating mass, hence constituent 
equations usually contain time derivatives of those tensors, and depend on a small 
number of viscosity parameters, to be added to the elastic Lame coefficients A and fi 
— see below. There is however no unique way in which this idea can be carried over 
from a one-dimensional system into a constituent equations set for a three-dimensional 
solid, and so different alternatives result in different models for the purpose. 

In this paper we propose to study and discuss the results of applying to a spherical 
GW detector the equations of various such phenomenological models, in view to 
determine how quality factors change from mode to mode in each case. This, it 
is hoped, will generate some insight into the nature of the viscous processes which 
take place in a specific spherical GW detector, and help to assess on the basis of 
spectral measurements which particular class of viscoelastic solid a given material 
belongs to. In turn, better understanding of material's macroscopic properties should 
also contribute relevant information to the currently important issue of spherical GW 
detector design. In section ^ we present the general equations, and in the subsequent 
sections we successively consider the Kelvin- Voigt, Maxwell, Standard Linear and 
Genaralised Mechanical models. As we shall see, the sphere's vibration eigenmodes 
always group into the usual families of toroidal (purely torsional) and spheroidal 
modes, but different quality factor dependences arise in different models. In section ^ 
we match our theoretical results to a number of experimental measurements of quality 
factors performed on various spheres, and discuss which is the most appropriate model 
to fit such data. In the final section we present a summary of conclusions, and two 
appendices are added to clarify a few mathematical technicalities. 

2. The general equations 



GWs bathing the Earth are known to be extremely weak — see e.g. |25 for a review — 
so that the classical equations of linear Elasticity are very good to describe the GW 
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induced motions of a spherical antenna in the expected frequency range, roughly 10^- 
10'' Hz. These equations are 

where s(x, t) is the field of displacements in the elastic body, and f(x, i) is the GW 
induced density of forces acting on the solid — see |l6| for full technical details. Uij is 
the the stress tensor, and is related to the strain tensor 



1 / dsi dsj 



2 \dxj dxi 



(3) 



through a set of constituent equations. In the case of a non-dissipative (ideal) solid 
these are simply the expression of Hooke's law. In a dissipative one, such equations 
include time derivatives of both and aij to account for internal friction effects. 
Constituent equations are of the following general type: 

L{sij , Sij , Sij , . . . ; CTij , a^j , CTy ,...) = (4) 

In this paper we shall limit ourselves to linear constituent equations, an excellent 
approach for a GW detector, as already stressed. In the simplest instance, only first 
order derivatives will appear in (^, and we shall consider this first. Then we shall 
also devote some attention to more complicted models. 

The equations of motion must be supplemented with suitable boundary 
conditions. We shall prescribe the usual ones 

aij Hj = at r = R (5) 

where n is a unit outward-pointing vector, expressing that the surface of the sphere 
is free from any tensions and/or tractions. 



3. Kelvin Voigt model 

This model assumes that the solid is homogeneous and isotropic, and is characterised 
by the following constituent equations p^ : 

A + A'^^ skk 5,j + 2 (^^ + ^'^^ s„- (6) 

The constants A and are the usual Lame coefficients describing the purely elastic 
behaviour of the body while the positive coefficients A' and /i' parametrise its 
viscous properties!]], which are proportional to the change rate of the strain tensor, 
dtSij. 

If equations are replaced into we obtain the equations of motion for s(x, t): 



(A + ^) + (A'+A.')| 



V(V-s) + f(x,t) (7) 



The solution to this system of coupled equations can be expressed in terms of a 
Green function, the construction of which requires explicit knowledge of its eigenmode 
solutions. As is well known, the latter correspond to the free oscillations of the solid 

II These coefficients are actually analogous to those which describe the viscosity of fluids in 
Hydrodynamics: shear viscosity (/i') and bulk viscosity (2/i' -|-3A'). 
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— no density of external forces in the rhs of (^. As usual, we attempt to find such 
eigen-solutions in the factorised form 

s(x,t)=T(i)s(x) (8) 

which results in the following (dots on symbols stand for time derivatives): 



pf{t) s(x) = ^T(t) + ^'T(t) V2s(x) 



i\ + fi)T{t) + {X' + fi')f{t) V(V-s(x)) 



(9) 



Next we use the well known decomposition of a three dimensional vector field into 
its irrotational a divergence free components pq] : 



s(x) = S;(x) + St(x) 



V-St(x) = V X s/(x) = 



(10) 



By the methods described in Appendix A, it can be seen that the equations 
satisfied by St(x), s;(x), and T{t) are (A.5)-(A.S), i.e., 



V^st + /C^ St = 



and 



(11a) 
(11&) 



= 







(12a) 
(12&) 



Since T{t) must fulfil both equations (116) and (12i), the separation constants K.^ 
and are not independent. The binding relationship is established after it is realised 
that the solution to those equations is of the form 

=7t 



where 7 is in general a complex quantity. We readily obtain 



(13) 



[A + 2^ + 7 (A' + 2^')] ' (14a) 
/C2= -p7'[M + 7/^r' (14&) 
The values 7 can possibly take on are determined by the boundary conditions, 
equations The reader is referred to Appendix B for a detailed description of the 
eigenvalue algebra of this problem. Just as in the non-dissipative case, there are seen 
to be two families of eigenmodes: toroidal (purely torsional) and spheroidal. Viscous 
effetcs are however small in practice, as inferred from the narrow linewidth of the 
measured resonances. This means that the following inequalities hold in any cases of 
interest to us: 



A' 1 
— < - 



(15) 



where oj is the frequency of the mode considered. We shall use these inequalities to 
estimate the roots of the eigenvalue equation ( B.15| ) (see Appendix B) perturbatively 
from the non-dissipative ones, already known — see ||l6| for full details. Clearly thus, 
our procedure will be valid for the lower frequency modes. We proceed sequentially 
for the two families of eigenmodes. 
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3.1. Toroidal modes 

These correspond to the solutions to (see ( B.15| )) 



/3i(/Ci?) = (16) 

which is formally identical to the toroidal eigenvalue equation for a non-dissipative 
solid sphere |l^. If we call fcj; the toroidal wave numbers of the latter, we have 

^nl = knl ^ \l (17) 

and hence, by equation (|14t 




= + -'1^ (18) 



where the last approximation depends on the validity of the assumption (|15|). 
Expression (|l^) nicely shows how this Kelvin- Voigt model predicts exponentially 
damped eigenmode oscillations. If we recall that such damping is expediently described 
in terms of a quality factor Q — see (|l|) above — then we discover that the prediction 
of the model is that 

Q^. = ^KO" (19) 

i.e., the quality factor for toroidal modes is inversely proportional to the frequency of 
the mode. The amplitudes of these modes have the form 

s^,,(x,i)=si(x,t)e--*/«, (20) 

/i to 

where the subindex KV stands for 'Kelvin- Voigt', while E refers to the standard 
frictionless case, whose amplitudes are those given e.g. in reference |p^ . 

3.2. Spheroidal modes 



A second alternative to find a non-trivial solution of the linear system ( |B.9| ) is to 
impose the condition 

PifaR,^^:^] f33{ICR)-l{l + l)l3i{QR)l3i{ICR)^0 (21) 



This is characteristic of the spheroidal eigenmodes. By virtue of equations (14q) 



and ( 146 ), this relationship can be translated into a condition to be fufiUed by 7, and 
which depends on the ratios A//i, A'/A and fJ-' / as well as on the multipole index I. 
In this case, as we are not dealing with an eigenvalue problem of a selfadjoint operator, 
complex solutions to equation (|2|) are allowed — indeed, expected. An exact solution 
of that equation implies a separation of its real and imaginary parts, followed by 
numerical calculations which determine the angular frequency and quality factor of 
the quasinormal mode at hand. We are interested in materials with long decay times, 
so we shall set up a perturbative solution to equation (|2l|), using 

(22) 
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as the small perturbative parameter. In other words, we assume that the 
approximation (^5|) holds. Here, w stands for a generic spheroidal eigenfrequency of 
the non-dissipative solid. Obviously, the unperturbed solution, i.e., that corresponding 
to e = 0, is the elastic solid's solution, already discussed in reference Q. We thus 
introduce the perturbative expansion 

7 = 70 + 7ie + O(e^), 70 = -icj, (23) 

and (j2^), we obtain perturbative expansions for the 
parameters /C and Q which can be written as 

/C = A:o + fcie + 0(e2), Q ^ q„ + + 0(6^), (24) 

where 



Using equations (|14cD, ( 



fco 



qo = q = uj 



P 



X + 2fi 



qi = I 



P 



X + 2fi 



71 



h' + 2uj 
h + 2 2 



(25 a) 



(25 &) 



In the above equations, k and q are the parameters appearing in the elastic 
sphere's case, and we have introduced the dimensionless ratios 



-a, 

P 



P 



(26) 



which are both zero order quantities. We can now perform the perturbative expansion 
of the eigenvalue equation. In order to ease the resulting expressions, let us introduce 
the following notation: 



l{l + I) f3i{K:R) 

PiiQR) 

P3{ICR) 

A + 7A' 



(34 QR, 



P + IP' 



l{l + l)/3i(fci?) + l{l + l)l3[{kR)kiRe 

Bo + BikiRe 

(5x{qR) + (5[{qR)qxRe 

Co + CiqiRe 

PsikR) + l3'^{kR)kiRe 

Do + DikiRe 



Pi{qR,h) 



(3',{qR, \lll)q^R~-{h' ~h)jMR) 



= AQ + {AiqiR-iA\)e 



(27a) 
(276) 
(27c) 

(27rf) 



where a prime on a function denotes differentiation with respect to its first argument. 
We note that the functions Aq, . . . , (defined in every second line of the above set 
of equations) are real. With this notation, the zero-th order form of equation ( pT| ) is 



A(,Do - CoBo - 



(28) 



which is simply the condition that w be a spheroidal eigenvalue of the purely elastic 
case. On the other hand, the first order expansion of (|l]) yields 



{AoDi - CoBi)kiR + (AiDo - CiBo)qiR = iA[Do, 



(29) 
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whence, using the form of ki and qi, the value of 71 ensues. It can be written as 

li^~'^[f{kR,h,h')]-' (30) 
where we define the dimensionless function / by|^ 

f{kR,h,h') = 

2A[Do{kR)-^ - AoDi + CqBi - (AiDo - CiBq) {h' + 2){h + 2)-3/2 ^ ' 

We note that the first order correction obtained for 7 is real. Therefore, to this 
order of approximation, the frequencies of vibration remain unahered, and are the 
same as those obtained for the elastic solid. Moreover, ki and qi happen to be purely 
imaginary. Therefore, the modulus of the radial functions appearing in the spatial 
part of spheroidal quasinormal modes of vibration will also be the same as those of 
the ideally elastic solid, for the corrections to k and q will just introduce, to first order, 
a complex phase factor. 

Summing up, we have shown that while the spheroidal normal modes of vibration 
of an elastic solid are given by an expression of the form ]l6| 

s;g(x,t) =e*'^"'* [AJir)Yirn{e,ip)n-Br,i{r)inxLYim{0,ip)] (35) 

the spheroidal quasinormal modes s^y of a Kelvin- Voigt solid are obtained from the 
normal modes of the elastic solid according to the following 



(x,i) 



iuj^it — uj^it/Q 



e*>^iMA„z(r)yi™(0,(^)n 



- e*x^MB„z(0^nxLyi™((?,(^)] (36) 
the qualitity factor being given by 

Q^i^4^f{kruR,h,h') (37) 

where, it is recalled, / is given by (HJ) as a function of the mode and the coefficients 
characterising the viscoel astic p rope rti es o f the s olid. The real phases Xi 2{^) '^^^ be 
computed from equations (B^)-(B^), (25q), and (25?). Nevertheless, the explicit (and 
cumbersome) form of these phases is largely irrelevant, and whe shall omit its explicit 

^ The case in whichf takes its simplest form is that of monopole modes. We know (see jl^) that 
when 1 = equation (Kfl) is no longer valid, and must be replaced by 

A + 7A' ^ 



Using now the expansion (27c) and the fact that Aq 
eigenfrequencies, we obtain for the first order correction to Q, 

qi = iA[A-^ 



(31) 

for the uperturbed monopole 

(32) 



and therefore, using the notation of equation ( |30[ ) and the relationships (25c)-(25f), 

-1 

f{qR,h,h') 



h' 



2 A 



h- 



gRAi 



(33) 



□ 
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Figure 1. Plot of the function f{kR,h,h') — see ( |34| ) — for a Poisson ratio 
fT = 1/3 (which means h = 2), and for the first 18 (spheroidal) modes of the sphere's 
spectrum and a few values of the viscoelastic ratio h' . Note that f{kR, h, h') = 1 
i{h = h'. 



form here |g7| . They merely introduce a position dependent shift in the phase of the 
vibrations which is of order e, therefore not hkely to give rise to measurable effects. 
More interesting and physically relevant is the behaviour of the function / giving the 
precise dependence of the quality factor on frequency. First of all, it is easily seen that, 
for the special case h = h', f is equal to 1, and thus the quality factor is proportional 
to as was the case with toroidal modes. But when the aforementioned equality 
does not hold, numerical calculations are needed. Figure |l| shows that this function 
has a rather irregular dependence on frequency. We have represented / for the first 
18 eigenvalues of the spheroidal spectrum, where the inequality (|5|) safely holds. It is 
clear from equation ( ^4[ ) that holding fixed h and kR leaves us with a linear function 
of h', whose slope varies from root to root. 



4. MeLxwell model 



In this section we shall consider constituent equations given by the so called Maxwell 
model. Like the Kelvin- Voigt, these equations only involve first order time derivatives. 
As we shall see, quite different predictions will be obtained for the quality factor 
dependence on the mode frequencies. 

The Maxwell model is also isotropic and homogenous, and is characterised by the 
following constituent equations |Q : 

d 

dtOij + a akkSij + 13 aij = — (A Skk S^j + Sij) (38) 

Here, the constants A and fi are again the Lame coefficients describing the elastic 
behaviour of the body, while the constants a and (3 parametrise the effects due to 



Resonance widening. 



9 



1 + - ) CTy (x) = ( A - a—^-:^—^ ] Sfcfc(x)(5„- + 2/xSy(x) (41) 



internal friction. To construct factorised solutions, we must also factorise both stress 
and strainQ: 

a,j (x, t) = e^* ay (x) , s,(x, t) = e''* s^{x), (39) 

The constituent equation is thus written, after separation of variables and 
contraction of its free indices, as 

(7 + 3a + (3) ffjj (x) = 7 (2/i + 3A) s^^ (x) (40) 

and hence we have the following relationship between the spatial parts of stress 
and strain: 

— (7)1 X = A — a „ 

7/ V 7 + 3a + /5, 

Like in the Kelvin- Voigt model, we shall be mainly interested in the case of small 

internal friction, so we shall assume 

A' A « 1 (42) 
I 7 I I 7 I 

whence the following constituent relationship results: 

ay(x) = A (^1 - sfcfe(x),5,^. + 2^ (^1 - Q Sy(x) (43) 
where we have introduced a new constant, (5, defined by 

6^'-^a + P, (44) 

so that we can take as the parameters characterizing the Maxwell solid the set 
consisting of the Lame coefficients A, ^, and the parameters (3, 5 which describe 
internal friction. 

Let us compare equation (^ ) with that of the Kelvin- Voigt model — cf. (^ — 
once the separation of variables has been performed: 

CTy (x) = (A + 7A') Sfefc (x) 5ij + 2(/i + 7/i') s,j (x) (45) 

Comparing equations (^3|) and (^5|), we observe that the solution of the Maxwell 
model can be carried out, as regards the spatial part of s, following the same method 
used in the previous section for the Kelvin- Voigt model. In fact, we can directly take 
the expressions there derived, and make the substitutions 

^l' — > , A' — > -A(57"2 (46) 

which transform equation ( ^ ) into ( ^ ) . Thus the form of the solutions and boundary 
conditions for a Maxwell viscoelastic sphere are those of Appendix B| with constants 



/C and Q now given by the following functions of the parameter 7: 

2 = h/T^ \l+T^. i^^ + P)] (47a) 




IC = zJ^{j+^] (47&) 

where the approximation ( |42| ) has been taken into account, and h = X/fi. 

Thus the two families of quasinormal modes of vibration are also present in this 
model, and we describe them in the following subsections. 

+ Due to the equations of motion (0), if s(x,t) is assumed to be separable in the fashion (fel) then 
the strain tensor is separable too, and due to the constituent equation it is easily seen that tne only 
possible time dependence is of the form exp(7t). 
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4.I. Toroidal modes 

As already discussed, the allowed values for 7 are again those making the linear 



system (B^) compatible, and there are two alternative ways to accomplish this. The 
first possibility yields purely tangential (Ct — Ci = 0) vibrations satisfying once more 
the condition 



/3i(/Ci?)-0 =^ lC=J^uj'^i (48) 

V M 

being a toroidal eigenfrequency of the elastic sphere^ Using the relationship 



between 7 and K, for a Maxwell solid given by equations (47Q)-( |47q ), we obtain the 
allowed values 7J; as 

ill = -^^li - f (49) 

Again, toroidal quasinormal modes have two fundamental properties: they have 
the same set of eigenfrequencies as the elastic sphere (to first order in the parameters 
describing internal friction, (3 in this case), and also exactly the same spatial part (for 
all values of the viscosity parameters). The only difference between Kelvin- Voigt and 
Maxwell solids as regards toroidal modes appears in the dependence of the quality 
factor on cj: as equation ( ^9| ) shows, the quality factor in a linear Maxwell solid 
increases linearly with frequency. We can express all these properties by means of the 
following formulae: 

sI,(x,i)=si(x,t)e---*/Q, = ^ (50) 

relating Maxwell quasi-normal modes of vibration, s^^(x, t), to elastic normal modes, 
s^(x, t), for the toroidal family. 

4-2. Spheroidal modes 

In order to handle the spheroidal family, we shall resort again to the perturbative 
expansions already used in the Kelvin- Voigt case, and also in the toroidal family, 
just described. The Maxwell model trivially reduces to the perfect elastic case when 
[3 — 5 — 0, hence we can take as the perturbative parameter 

(51) 

where uj is the elastic eigenfrequency to which 7 approaches when both (3 and 5 



approach zero. Perturbative expansions in the fashion of section 3.2 can now be 
introduced: 

7 = -iw-|-7ie, IC = k + kie, Q — q + qi£ (52) 



where the first order corrections ki and qi are given by equations ( ^:7q )-( |47q ) as 
functions of 71 : 

fci = (71 + I) (53a) 

(536) 
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The zero-order ratio h' is now given by 

h' = ^h (54) 

With this definition, together with that of the perturbative parameter, the 
expressions at hand are formaUy identical to those of the Kelvin-Voigt model, 
and therefore the solutions to the Maxwell model share all their properties with 
their Kelvin-Voigt counterparts; the exception is the dependence of the quality 
factor on frequency: the product 71 e, which gives the exponential decay, is now 
independent of lu. 

Summing up, spheroidal quasinormal modes of the Maxwell solid, s£f(x, t), are 
related to spheroidal normal modes of a perfectly elastic sphere by the equations 

sU^,t) = sg(x,i)e--*/Q+>^W , Q^^ = ^f^kr,iR,h,h') (55) 

where the function f{kR, h, h') is again given by (^4|). 

As we see, the only difference between the behaviour of Maxwell and Kelvin-Voigt 
viscoelastic solids, when the internal friction effects can be considered small, appears 
in the dependence of Q on frequency. We must however stress that, under other 
conditions (e.g . static load), both models show larger divergences in their physical 
properties [E8|. 



5. Other Models 



In this section we review other models which have been proposed to address the 
dynamics of a viscoelastic solid. They are generalisations of those in the two previous 
sections. We shall however not attempt to find complete solutions to all of them, as 
it eventually becomes too cumbersome. We shall however discuss in this section some 
of their most relevant traits. 



5.1. The Standard Linear Model 

The Standard Linear Model (SLM) for a viscoelastic solid is a generalised combination 
of the Kelvin-Voigt and Maxwell models. The constituent equations take here the 
form: 



where the effects of internal friction are described in this case with the aid of four 
constant parameters: a, (3, a', and /?'. When looking for factorised solutions, the 
equations of motion and the above relationship force a time dependence of the form 
e'*'* for both stress and strain. When such a dependence is introduced in equation (p6|), 
we obtain the following relationship between the spatial parts of the stress and strain 
tensors: 



(l + 27/3)a,,(x) = 



A + a'7 — 



3A + 2/i + 7(3a' + 2/3)' 



l+7(3a + 2/3) 



Sfefc(x) Sij 



+ 2(/i + /3'7) 5y(x) 



(57) 
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The case of small internal friction is treated by first order approximation in the 
quantities parametrising viscous processes, i.e., 

a, /?, a', /3' < I 7 (58) 

When such approximation is made, the above equation reduces to 

cry (x) = (A + A'7)sfefe (x)5y + 2(/i + //7)s.y (x) (59) 

where we have introduced two new constants given by 

H' = P' ~ 2/3/i , y = a' - 2X13 - (3A + 2^i)a (60) 

Therefore when equation (|5^) holds, the SLM reduces to a Kelvin- Voigt model, 
i.e., for small internal friction, both models have the same set of quasinormal modes 
of vibration, which are characterized by the constants A, pt, A' and /i', the latter being 
given, for the Standard Linear solid, by equations (|60[). 



5.2. Generalised mechanical models 

The models analyzed so far are the simplest ones obtained by three dimensional 
generalisations of mechanical viscoelastic models composed of linear springs. They 
give rise to differential constituent relations, with time derivatives up to the first order. 
Considering more involved generalisations yields differential relations involving higher 
order time derivatives of strain and stress — see e.g. Thus, quite independently of 
any reference to the underlying mechanical model, we can consider general differential 
relations between stress and strain including any number of time derivatives. To ease 
the formulation of such differential constituent equations for the case of isotropic 
and homogenous solids, we shall introduce the trace-free parts of the strain and 
stress tensors, s^^, (usually termed deviatoric components in the literature on 
viscoelasticy |^^), and their traces, s and a {dilational components), defined by 

^'ij = ~ 3 S'^y" ; s = Skk (61a) 

'^'ij = ^ij - -^f^^y ' (T = CTfcfc. (61 &) 

In terms of the above quantities, the linear Hooke law for an elastic solid takes 
the form 

a=(3A + 2^)s <=2Ai4. (62) 

while the constituent equation of an SLM is written 



1 + (3a + 2/3) 



dt 



(3A + 2fi) + (3a' + 2(3') 



dt 



1 + 2/3 



dt 



ai, ^{2^l + 213' 



■'1 
dt 



(63 a) 



(63&) 



This equation can now be generalised to include higher order time derivatives. 
We can thus consider viscoelastic models whose constituent equation is given by 



R{dt)a = S{dt)s , i?'(9t)4 = S'{dt)s',, 



(64) 
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where we have introduced formal polynomials 

N-l N-1 

R{x) EE nx^ , R'{x) = ^ r'lX^ (65a) 

N-l N-l 

S{x) = six^ , S\x) = Y 4"^^ (65&) 

so that a general differential model is given for each set of 47V real constants rj, r[, 
si, and s'l characterising the solid. Some of these constants may vanish. The Kelvin- 
Voigt, Maxwell and SL models considered above are of course special cases within 
this general class. Several procedures have been proposed in the literature to solve 
the general equations — see pO| for a review and further reference — , which can be 
applied to the solid viscoelastic sphere problem. We now sketch how they work in this 
case of our interest. 

Let dij (x, il) and Sij (x, fi) be the Fourier transforms of the stress and strain 
tensors, and Si(x, fi) that of the displacement vector field: 

/oc 
cry(x,t)e^"*dt (66a) 
-OO 

/OO 
s^{x,t)e'"Ut (666) 
-OO 

/OO 
s^j{x,t)e'^Ut (66 c) 

-oo 



In terms of these, Fourier transforms the constituent equations (|6J) read 
Riifl) , ^, _ R'iifl) , 

while the equations of motion are 



^^_l(R{in) R'{in)\ R'(i^) 



Comparing the above equations with the corresponding ones for normal modes 
of vibration of elastic solids, and the constituent relation ( |67| ) with |6^, we note that 
the problem of finding solutions to the equation of motion of a general viscoelastic 
differential model reduces to that of finding the normal modes of vibration of an elastic 
solid having complex Lame coefficients given by 

~ 1 (Rm R'm \ 1 R'm 

= 3 \-m -^)) ^ ^ 2^ 

where the allowed values of are obtained as the solutions to the elastic solid's 
eigenfrequency equation when the above complex coefficients are used instead of the 
real, constant Lame coefficients A, /i. Generally, 17 will have complex values, thus 
giving rise to damped system oscillations. After solving for fi, the spatial part of 
the solutions is obtained from that of the normal modes by simply substituting the 
old, real-valued constants w, A and by the new complex values fi, A and jl. This 
method for solving the viscoelasticity is often termed in the literature on the subject 
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the Correspondence Principle [ pO[ , and as a matter of fact our previous derivations of 
the form of the quasinormal modes for Kelvin- Voigt, Maxwell and SL models can be 
seen to be special cases of its application. The method is applicable to any boundary 
value problem whose elastic counterpart is solvable. The case of small internal friction 



(i.e., first order approximation in the coefficients of the polynomials (S5q)-(55?) has 
been considered by Graffi | |30| for one dimensional wave propagation. 

The three dimensional spherical case is also solvable, as we know. The toroidal 
modes are relatively straightforward to obtain from their elastic counterparts due to 
the simple form of their eigenvalue equation, while the spheroidal ones demand more 
complex algebra, which becomes increasingly cumbersome as the order N of the model 
increases. We shall present here the general solution for the toroidal modes for any 
differential viscoelastic model, whereby we shall obtain the dependence of their Q on 
frequency. This will also be the approximate dependence for the spheroidal modes, if 
friction effects are small, as was the case with the first order models analysed so far. 
A complete solution for the latter modes can also be systematically found, but will be 
omitted due to its scarcely useful algebraic complexity ]2^ . 

5.2.1. Toroidal modes As discussed above, the boundary equation for the toroidal 
modes in a general viscoelastic model is obtained from the eigenvalue equation of the 
elastic model: 



/3i(fci?) = 0, k^J-uj (70) 
V 

Upon substitution of fi by /i, we obtain 



We know that the only solutions to the eigenvalue equation [t^ are the real 
eigenfrequencies of the elastic sphere uj'^i, and therefore the allowed values for are 
given by the implicit relationship 



V"rW"-""' ^^^^ 

Let us now write the polynomials S' and R' in the form 

N-l / N-\ \ 

S\x) = 1 + e ^ s\x' , R;{x) - 2/i 1 + e ^ r[x' (73) 
;=i V /=i / 

as a suitable one to imply that internal friction effects are small, letting 

e < 1 (74) 

The quantities rjw' and sjw' are thus zero order in e and dimensionless, to being 
a toroidal eigenvalue of the elastic case. We then introduce an expansion for 
in the small parameter e, whose zeroth order term corresponds to a given toroidal 
eigenfrequency lo of the elastic solid: 

rj = u; + el7i (75) 
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Under the above conditions, we have 

+ «,^.-'(.l-l)/2 (76) 

and the value of fli foUows when introducing the above expansion into equation 
yielding 

1=1 

In terms of Qi, the quality factor reads 

^ e Ini[51i] ^ ^ 

where Iin[-] denotes the imaginary part of its argument. Thus we observe that, 
as regards toroidal modes, using a general differential model gives us a polynomial 
in w for l/Q, with no independent term, so that constant Q is not allowed by these 
models. The polynomial only contains odd powers of the unperturbed frequency uj. 
In general, whenever ti^O for even /, the real part of fii will not vanish, and the 
angular frequency of the periodic component of the quasinormal modes shall undergo 
first order corrections. Hence, in order to preserve the elastic spectrum to first order, 
our model must satisfy the conditions 

ti^O {I even) (79) 

Provided the preceeding equation holds, the corrections to K- will be purely 
imaginary, and therfore the modulus of the spatial part of the modes will remain 
unaltered, the only effect of viscosity being the addition of a point dependent phase 
in the fashion of equation (p5|)p|. 

The calculation for the spheroidal quasinormal modes can be performed along the 
same lines but, as we have seen, the algebra is considerably more involved already in 
the simplest models. It does naturally become more cumbersome as the order of the 
model increases, so we omit a detailed discussion of its technicalities here. 



6. Experimental data 

We now present a confrontation of the above theoretical models with available 
experimental measurements, some of them published [^^, and others unpublished 
as of this date |^ . 

It is very important, in order to rightly assess the suitability of the theoretical 
models, to keep in mind that measurements of Q are intrinsically difficult. The most 
important source of problems is the fact that almost any experimental environment 
conditions (suspensions, readout, electronics, etc.) significantly couple to the primary 
oscillator (the sphere itself), and almost invariably result in a degradation of the 
measured Q. The way and intensity in which this happens is however utterly 
unpredictable. As a rule, theoretical Q values should thus be expected larger than 
measured ones. 

* This correction will only appear provided that N >3. This is why it was absent in the toroidal 
families of the previously discussed models. 
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Figure 2. Quality factors of a sample sphere of A15056 (top), and a 7% alloy 
of CuAl (bottom). The first had a diameter of 153 mm, and measurements were 

taken at 80 mK, while the second was 150 mm in diameter and had a temperature 
of 20 mK. As can be appreciated, the Q's of the lower frequency modes match 
the theoretical predictions fairly well, assuming a Maxwell model with viscoelastic 
ratios h' = I and h' = 0.08, respectively. The higher mode in CuAl, however, 
considerably deviates from calculations. 



Then, the mode frequencies slightly split up as an unavoidable consequence of 
the necessity to suspend the sample in the laboratory a fact not considered in the 
above models. This, at times, generates a forest of frciquency peaks when different 
multipole modes happen to be close to one another. For example, within a mere 0.8% 
of their nominal value, we find the three fundamental frequencies: wf^, and 
which, if split up, result in as many as 23 Fourier peaks in a very reduced frequency 
interval. On the other hand, however, the Q's of the various multiplet components 




Figure 3. Quality factors of a sample sphere of steel; the sample had a diameter 
of 107 mm and measurements were taken at a tempreature of 4 K. Lower frequency 
Q's are fairly well fit by a Maxwell model with a viscoelastic ratio h' = 0.3, while 
higher frequency Q's require a Kelvin-Voigt with viscoelastic ratio h' = 2.5. In 
this case there more experimental points than theoretical due to mode splitting 
caused by the sphere suspension — see text. 



often group around clearly distinct values, close to those predicted by theory. In 
such cases, therefore, the viscoelastic model is helpful to tell the different multiplet 
members from one another, an otherwise very difficult task. 

Rather than characterising a given viscoelastic model by the two parameters A' 
and or a and f3 — see sections || and || — , it is expedient for the purposes of 
this section to use a different parameter pair. This will be chosen as the Q of the 



lowe st f requency toroidal mode, Qo = (5i2, and the viscoelastic ratio h' of sections 3.2 



and |4.2| . The following expressions thus hold: 



Kelvin-Voigt: { ) ' (80a) 

f{kPR,h,h')Qo 



Maxwell: { ^ (806) 



obviously a rewrite of the corresponding ones in the respective sections. 

In figure we plot the measured versus theoretically calculated values of Q for 
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two sample spheres of Al 5056 and a 7 % alloy of CuAl — see more technical data in 
the caption to the figure. Agreement is quite good, when a Maxwell model is used, 
except for the higher frequency mode, reported in the case of CuAl. It appears that 
a better fit can be accomplished when a combination of both Maxwell and Kelvin- 
Voigt models is used, such as we see in figure ^, where Q's of a sample sphere of steel 
are considered. The latter plot contains more experimental points than theoretical, 
and we see here a clear example of frequency splittings due to symmetry breaking 
caused by suspension, normally a diametral or semi-diametral bore drilled across the 
sphere [ plf . 

7. Conclusions 

In this paper we have addressed the problem of whether it is possible to systematically 
characterise the linewidths or, equivalently, the Q's of the oscillation eigenmodes of a 
given spherical GW detector. To this end we have considered various phenomenological 
models, selected from the specialised literature on the subject, and solved the equations 
of motion in the cases of our interest. Different models are seen to predict different 
frequency dependences of the quality factors for the lower modes, which are the ones we 
have paid attention to, and the ones relevant for GW detection purposes. For example, 
in a Kelvin- Voigt solid the Q of a given mode appears to be inversely proportional to 
its frequency, while in a Maxwell solid it is directly proportional to it. 

As we have seen, however, the behaviour of a given material is generally not well 
described by a single such model but rather by a suitable combination of e.g. Maxwell 
and Kelvin- Voigt models, depending on the range of frequencies considered. This 
is not too surprising, as these models are phenomenological, i.e., they are not based 
on a detailed analysis of the physics of the various dissipation mechanisms which 
are responsible for the loss processes |3^, but on more or less plausible (though not 
unique) generalisations of the "friction is proportional to velocity" principle, which so 
accurately holds in a simple unidimensional harmonic oscillator. 

In spite of these limitations, and also in spite of the experimental difficulties 
inherent in the determination of Q, it appears that our treatment of the problem 
gives results which can be considered quite good, in view of the above limitations. A 
significant bonus of the present analysis is the possibility it offers to identify nearby 
frequencies in the sphere spectrum by their distinct value of Q, as predicted by theory. 
This is useful in GW detection science since, for example, the second spheroidal 
quadrupole frequency, w^^, is very important due to large cross section for GW energy 
absorption in that mode |^ yet it is very close (in frequency) to other modes with 
different Q's. 

Improvements on the results presented in this paper are certainly possible. But 
we feel they would have to be based on a different methodology, with more emphasis 
ab initio on Solid State lore, and/or on Materials Science. 
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Appendix A. 

A well-behaved three-dimensional vector field s(x) can be expressed as the sum of an 
irrotational, s/(x), and a divergence free, St(x), vector fields, respectively called the 
longitudinal and transverse components of s(x) p6[ |: 

s(x) = s,(x) +st(x) , V-Sf==Vxs,=0 (A.l) 

We now replace this decomposition into equation (|^) to find 



pfit){st + si) ^ liT{t) + p.'t{t) V\st+si 



(A + /i)T(t) + (A' + ^')r(t) V(V-sO 



Taking the rotational of this equation. 



V X 



= 



(A.2) 



(A.3) 



The vector between square brackets is thus divergence-free and irrotational, so it 
vanishes. We have therefore: 



St 



(A.4) 



Since the left hand side of the above equation does not depend on time, the term 
between braces in the right hand side must equal a (complex) constant, say —K?. 
Thus, 



iT + fi' f + IC^'^f = 



(A.5) 
(A.6) 



An analogous procedure, after taking the divergence of equation (A.2), gives us 
the corresponding formulae for the longitudinal part: 



V^s, + Q2 s; = (A.7) 
(X + 2fi)T+{X' + 2n')f+Q-'^f^0 (A.8) 

where stands for another complex separation constant. 
Appendix B. 

We describe in this appendix the algebraic operations which lead to the solution to the 
eigenvalue problem in a viscoelastic sphere. Equations (||) ought to be solved, subject 
to the boundary conditions (^. The latter can be cast in explicit vector form: 

(A + 7A') [V-s(x)] n + 2{fi + -ffi') (n- V)s(x) + 2{fi + jp')nx [V x s(x)] = (B.l) 

The irrotational and divergence free components, St(x) and s;(x), can be 
expressed by means of auxiliary functions ^!>(x) and V'(x): 
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Si(x) = Q-iCoV0(x) 

st(x) =i/C"iCiVxL7^(x)+iC2L?/'(x) 



(B.2) 
(B.3) 



where Cq, Ci and C2 are (so far) undetermined integration constants, and L = — ixx V 
is the "angular momentum" operator. Upon substitution of (B.2)-(B.3) into (0) it is 



readily seen that the functions (j) and ip are themselves also solutions to corresponding 
Helmholtz equations, 



(x) + Q^0(x) = , and + /CV(x) ^ 



(B.4) 



They have therefore the general form, using spherical coordinates {r,9,ip) for the 
vector x, 

0(x) = jiiQr) Yi„,{e, if) , and V(x) = ji{JCr) ip) (B.5) 

where ji are spherica l Be ssel functions of the first kind and Yim are spherical 
harmonics. The solutions (B.5) are those possessing regularity properties in the whole 
interior and boundary of the solid. We thus have 



V(l) 

VxL-0 



dr r 
-l{l + l)^-^^YUe,^)n 



r dr 
: jiilCr) iLYim{0,ip) 



in X LYi„i{e,ip) 



(B.6) 

(B.7) 
(B.8) 



Thes e expressions ought to be substituted now into ( |B.2| )-( |B.3| ), and then 
into (BJ) — recall that s = Sf + s;. It is found thct these are equivalent to the following 
homogeneous linear system 



/ 

V 

where 



-PliQR) 




-;(? + l)|/3i (/Ci?) 
§/33(/Ci?) 










Jiiz) 



d 

dz 
d^ 
dzJ 



Jliz) 



[Ji{z)] 



i/32(z) + {^-l} Mz) 



Pi{z,A)=l32{z) 



A 



(B.IO) 
(B.ll) 

(B.12) 
(B.13) 
(B.14) 
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The system (B.9) is to be satisfied by the constants Cq, Ci and C2, but has 
no meaningful solution unless the system matrix is singular, i.e., if its determinant 
vanishes. Therefore 



/34 (Qi?,^) + l)|/3i (/Ci?) 



Pi(ICR)det \ ^n-" 'A'+TA'V ^ 'Q^'^ ' =0 (B.15) 

V -/3i(s^) ^um ) 

This is an equation for the parameter 7, on which K, and Q depend through 



equations ( |l4q ) and (146). Clearly, there are two families of solutions, or eigenmodes, 
to (B.15) associated to the vanishing of either of the two factors in its Ihs, i.e., 
Pi{JCR) or the determinant of the displayed 2x2 matrix. They are called toroidal 
and spheroidal solutions, respectively. 
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